UNCLASSIFIED 


AD  NUMBER 

AD815126 

NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release,  distribution 
unlimited 


FROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  and  their  contractors; 
Administrative/Operational  use;  Mar  1967. 
Other  requests  shall  be  referred  to  Air 
Force  Systems  Command,  Space  Systems  Div, 
Los  Angeles  CA  90009-2960. 

AUTHORITY 

DoDD  5230.24,  18  Mar  1987 


THIS  PAGE  IS  UNCLASSIFIED 


Best  Available  Copy 


Aerospace  Corporation  No 
TR-i001(2240-30)-9 


Air  Fwrce  Report  No. 
SSD  TR-67-64 


THE  EFFECT  OF  RISE  TIME  ON  CRITICAL  DYNAMIC  LOAD 
FOR  A  SHALLOW  ARCH 

Prepared  by- 

Malcolm  H.  Lock 
Solid  Mechanics  Department 
Applied  Mechanics  Division 


Laboratory  Operations 
AEROSPACE  CORPORATION 


March  1967 


Prepared  for 

BALLISTIC  SYSTEMS  AND  SPACE  SYSTEMS  DIVISION 
AIR  FORCE  SYSTEMS  COMMAND 
LOS  ANGELES  AIR  FORCE  STATION 
Los  Angeles,  California 


FOREWORD 


4 


This  raport  is  published  by  tha  Aerospace  Corporation,  El  Stgundo, 
California,  under  Air  Force  Contract  No.  AF  04(699) -1001. 

This  raport,  which  documents  research  carried  out  for  Aerodynamics 
and  Propulsion  Research  Laboratory  from  1  August  1966  through 
1  December  1966,  wee  submitted  on  29  March  1967  to  Gapt.  John  T.  Allton, 
S8TRT,  for  review  and  approval. 


Approved 


Erector 
Engineering /Sciences  Subdivision 
Applied  Mechanics  Division 


Aerodynamics  and  Propulsion 
Research  Laboratory 


Publication  of  this  report  does  not  constitute  Air  Force  approval  of 
the  report's  findings  or  conclusions.  It  is  published  only  for  the  exchange 
and  stimulation  of  ideas. 


-ii 


ABSTRACT 


The  critical  ramp  loads  required  to  produce  snapping  of  a  shallow 
sinusoidal  arch  are  investigated.  Calculated  results  are  presented  for 
two  specific  arch  geometries.  These  results  illustrate  the  influence  of 
the  load  rise  time  on  the  level  of  the  critical  load. 
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*  generalised  coordinate  denoting  amplitude  of  m^1 
component  of  the  nondime  ns  ional  displacement  T|(?,  t) 

a  elastic  modulus  of  the  arch  material 

*  nondimenaional  geometric  parameter;  see  Eq.  (2). 

*  second  moment  of  area  of  the  arch  cross-section 
a  thickness  of  the  arch  cross-section 

*  radius  of  gyration  of  the  arch  cross-section 
>  span  of  the  arch 

a  pressure  load  acting  upon  the  arch 
■  (p/EIk)  (L/tt)*,  nondimens  ional  pressure  load 

*  critical  value  of  nondimensional  pressure  load 

*  /e ,  critical  static  buckling  pressure 

a  parameters  in  psramstric  resonance  analysis;  see 
Appendix 

a  time 

a  transverse  displacement  of  middle  surface  of  the  arch 
meaaured  from  the  baseline  of  the  arch 

a  Cartesian  coordinate  along  span  of  the  arch 

a  (i/k)  w  nondimensional  displacement  of  the  arch 

a  ttx/L,  nondimensional  Cartesian  coordinate 

«  mass  density  of  the  arch  material 

»  tfrr/L)*  (EI/p#h)*^,  nondimensional  time 


NOMENCLATURE  (Continued) 

nondimensional  rise  time  of  the  ramp  load 

2tt/< *>j,  nondimensional  period  of  free  vibration  of 
the  fundamental  symmetric  mode 

2  1/2 

(1+e  /2)  ,  nondimensional  frequency  of  free  vibration 

of  the  fundamental  symmetric  mode 


SECTION  I.  INTRODUCTION 


The  critical  step  pressure  loads  required  to  produce  snapping  of 
a  shallow  sinusoidal  arch  were  previously  determined  by  the  author  for 

a  wide  range  of  arch  geometries  .  These  criti  il  pressures,  which 
were  in  some  cases  as  low  as  73  percent  of  the  corresponding  critical 
static  buckling  pressures,  result  from  a  particularly  severe  dynamic 
load  condition  since  the  full  level  of  the  loading  is  applied  instantaneously. 
It  is  anticipated  that  the  severity  of  the  loading  will  be  reduced  in  those 
cases  wherein  a  finite  interval  of  time  is  required  for  the  load  to 
develop  to  its  maximum  level.  The  present  note  describes  an  investi¬ 
gation  of  this  effect. 


Lock,  M.  H.  'Snapping  of  a  Shallow  Sinusoidal  Arch  Under  a  Step 
Pressure  I  *d,'  AIAA  J.  4(7),  1249-1256  (July  1966). 
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SECTION  U.  ANALYSIS 

We  consider  a  shallow  sinusoidal  arch  subject  to  an  external  pres¬ 
sure  load.  The  temporal  behavior  of  the  loading  is  ramp-like  (see 
Figure  1)  and  the  spatial  distribution  is  sinusoidal.  The  load  is  charac¬ 
terized  by  a  rise  time  Tf  and  by  a  maximum  level  q^.  The  nondimen- 
sionalized  pressure  q(g*T)  is  written 

q(5.f)  -  sin  5  (1) 

The  nondime  ns  ional  transverse  displacement  T|(5,t)  of  the  arch  is 
approximated  by  the  following  modal  expansion 

T]{?.t)  =  e  sin  5  +  a^r)  sin  %  +  a^-r)  sin  2£  (2) 

where  the  generalized  coordinates  a^(T)  and  a2(T)  denote  the  amplitude  of 
the  first  symmetric  and  antisymmetric  modes  of  free  vibration;  the 
term  e  sin  g  describes  the  initial  shape  of  the  arch.  The  geometric 
parameter  e  is  defined  by  e  =  initial  height  of  the  arch  /radius  of  gyration 
of  the  arch  cross-section.  The  generalized  coordinates  a^f-r),  a are 
governed  by  a  pair  of  coupled  nonlinear  ordinary  differential  equations. 
Except  for  the  change  in  the  time  history  of  the  load  3j(t)»  these  equations 
are  identical  to  those  treated  in  Reference  1. 

The  relief  in  the  severity  of  the  loading  (as  compared  to  the  step 

pressure  case)  is  exhibited  by  displaying  the  variation  of  the  critical 

* 

load  q^  with  the  rise  time  parameter  a  where 

«  =  (1  +  Tr/‘rj)"1 

and  Tj  denotes  the  period  of  the  fundamental  symmetric  mode^. 

1 - - - 

The  rise  time  parameter  or  varies  between  the  values  of  zero  (static  load 

limit  r  /t4  -ao)  and  unity  (step  load  limit  t  /t.  -*0). 
r  1  r  l 


UMM.M 


ft 

The  critical  loads  were  determined  by  numerical  integration  of 
the  nonlinear  equations  of  motion.  These  loads,  which  are  the  lowest 
loads  at  which  the  response  of  the  arch  encompasses  the  corresponding 
snapped  equilibrium  state  of  the  system  (Reference  i)  ai  e  revealed  by 
a  "jump*  in:  the  level  of  the  response.  Results  were  obtained  for  two 
specific  arch  geometries;  namely,  e  =  4  and  e  =  7.  These  geometries 
were  selected  upon  the  basis  of  the  results  obtained  in  Reference  1  and 
are  representative  of  the  two  regimes  of  snapping  mechanism  that  were 
reve'led  in  that  study.  The  lower  value  of  e  corresponds  to  the  "direct* 
snapping  regime;  the  higher  alue  corresponds  to  the  "indirect"  snap¬ 
ping  regime. 

The  calculated  results  are  shown  in  Figure  2  in  terms  of  the  critical 
load  ratios  q^/q^.  where  denotes  the  corresponding  static  buckling 
load.  These  ratios  are  plotted  against  the  rise  time  parameter  er.  The 
results  obtained  for  the  two  geometries  are  seen  to  be  markedly  different. 
The  critical  load  ratio  for  the  case  e  =  4  increases  steadily  as  a  is 
decreased  until  it  reaches  a  maximum  value  of  unity  when  a  becomes 
zero.  Thus,  for  e  =  4,  the  level  of  the  critical  ramp  load  tends  mon o- 
tonically  to  the  corresponding  static  buckling  lead  as  the  rise  time 
ratio  t^/tj  is  increased.  In  the  case  of  e  =  7,  the  variation  of  the 
critical  load  ratio  with  a  is  marked  by  the  appearance  of  a  "jump*  in 
the  value  of  the  load  ratio  for  or  »  0.  57  (i.  e.,  for  Tr  /r^  m  0.  75)  and  by 
the  appearance  of  load  ratios  in  excess  of  unity  for  a  <  0.57. 
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Figure  2.  Variation  of  Critical  Load  Ratio  With  Rise  Time. 


SECTION  in.  DISCUSSION 


The  different  character  of  the  critical  load  results  may  be  attrib¬ 
uted  to  the  different  mechanisms  of  snapping  that  are  operative  in  the 
two  cases.  As  noted  previously  the  critical  loads  for  e  =  4  are  con¬ 
trolled  by  the  "direct"  snapping  mechanism.  The  arch  response  is 
primarily  in  the  symmetric  mode  (Figure  3)  and  snapping  occurs  when 
the  amplitude  of  the  generalised  coordinate  a^(r)  coincides  with  an 
unstable  equilibrium  state  (Figure  4).  For  e  =  7,  the  occurrence  of 
snapping  is  controlled  by  a  complicated  interaction  between  the  sym¬ 
metric  and  antisymmetric  modes  of  the  system.  A  more  detailed 
discussion  of  the  two  mechanisms  of  snapping  will  be  found  in  Reference  1. 
Briefly,  the  antisymmetric  mode  is  parametrically  excited  by  the  initial 
response  in  the  symmetric  mode.  If  this  excitation  is  sufficiently 
strong  the  interaction  from  the  antisymmetric  mode  back  to  symmetric 
mode  will  precipitate  snapping.  This  mechanism  of  snapping  gives 
rise  to  the  type  of  supercritical  response  shov.v*  In  Figure  5.  The 
presence  of  a  "jump”  in  the  critical  load  ratios  and  the  appearance  of 
load  ratios  in  excess  of  unity  may  be  associated  with  the  role  that 
parametric  resonance  plays  in  this  interaction.  A  simple  analysis  to 
illustrate  the  role  of  parametric  resonance  is  presented  in  the  Appendix. 

In  conclusion,  the  presented  results  give  an  indication  of  the  relief 
in  critical  load  level  that  is  introduced  by  a  finite  load  rise  time.  The 
calculations  show  that  the  critical  load  levels  required  to  produce 
snapping  are  increased  from  the  corresponding  step  load  values  for 
geometries  wher*  dynamic  weakening  occurs  (i.e.,  critical  step  loads  < 
critical  static  leads).  The  increase  is  gradual  if  the  snapping  process 
is  controlled  by  the  adirect"  snapping  mechanism,  the  increase  is  , 
more  rapid  if  the  "indirect"  snapping  mechanism  is  operative.  In  the 
latter  case  the  dynamic  weakening  effects  have  disappeared  for 
t?/tj  <  0.  75. 


- -  STATIC  LOAD  DEFLECTION  CURVE 

-  LOCUS  OF  MAX.  DISPLACEMENT 

UNDER  RAMP  LOAD  WITH 
Tf/r|  =  0.3 


Best  Available  Copy 


Figure  5.  Response  of  Arch  to  Supercritical  Ramp  Load: 
Zero  Damping  and  e  ~  7. 


Figure  5.  Response  of  Arch  to  Supercritical  Ramp  Load: 
Zero  Damping  and  e  *  7. 


APPENDIX 


PARAMETRIC  RESONANCE  ANALYSIS 

The  role  of  parametric  resonance  in  the  snapping  phenomenon 
may  be  illustrated  by  the  following  analysis. 

If  the  interaction  of  the  antisymmetric  mode  response  is  neglected, 
an  approximate  solution  for  the  generalised  coordinate  a ^  (t)  is 

(3) 

(4) 

(5) 

(6) 

and  Where  0  denotes  a  constant  phase. 

The  frequency  is  given  approximately  by  the  relation  (see  Refer¬ 
ence  1  for  method  employed) 


where 


al  **  —j—  1 1  -  6  cos  (fa>T  +  0)  J 

■j-M) 


6  =  | sin 


T  „  2n 

Ti  -  srr 
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If  the  solution  (3)  '.§  substituted  Into  the  linearised  version  of  the  equation 
governing  a2(*)  we  obtain 

,2 

4  a2  * 

~  yr  +■  (R  +  S  cos  f)  a,  B  0  (8) 

dtz  3 

where 

t  *  (wr  +  0)  (9) 

R-(=)2[“f]  <“> 

S  *  *(s)l(^-)  I1') 

Eq.  (8)  may  be  recognised  a«  the  Mathieu  equation;  two  of  the  sones 
of  parametric  resonance  determined  by  th)  s  equation  are  shown  in 
Figure  8  (sones  A  and  B  In  the  figure).  S-ich  sones  define  regions  where 
the  solutions  of  Eq*  (8)  exhibit  an  exponential  growth  of  the  form  e^T. 

The  iso-p  lines  shown  In  the  figure  conne  ct  points  of  equal  growth  rate; 
the  boundaries  of  the  sons  corresponding  to  p  =  Q.  Considering  these 
lines  as  a  measure  of  the  strength  of  the  parametric  excitation*  we  see 
that  the  excitation  is  stronger  In  the  interior  regions  of  the  sones.  The 
figure  also  shows  straight  line  loci  of  the  values  of  the  load  ratio  q  ^  Ajj 
for  fixed  values  of  t.,/T|*  These  loci  intersect  the  senes  of  parametric 
resonance.  As  q^/^j  ta  increased  from  sero,  for  fixed  the  loci 

first  intersect  sons  A.  If  the  excitation  in  the  region  of  the  rone  traversed 
by  the  loci  is  sufficiently  strong >  snapping  will  occur  at  some  critical 
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gure  6.  Zones 


value  of  the  loadratio.  However,  as  /tj  is  raised  the  regions  ofthe 
soue  tirayersefliby  lppi  become  increasingly  weaker,  ip  excitation, 
Eventually  the  excitation  ip  this  none  becomes  top  weafc  to  precipitate 
■napping.  A#  a  result,  the  load, ratios  must  be  increased  above  unity  to 
•atlPly  the  condition* erf  parametric  resonance  in  aone  B,  the  snapping 
phenomenon  now  being  precipitated  hy  the  excitation  in  this  zone.  Thus . 
It  la  seen  thut.  the  change  Ip  the  aone  ad, parametric  resonance  accounts 
<«?  both  the  jump  ln,critical  loed  ratln  and  the  appearance  o*  load. ratios 
greater  than  unity,  (figure  ?)« 
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